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Abstract 

We deduce approximate equations of state from experimental measurements 
in pp and pp collisions. Thermodynamic quantities are estimated combining 
the measure of average transverse momentum (p?) vs pseudorapidity density 
with the estimation of the interaction region size from measures of Bose 
Einstein correlation, or from a theoretical model which relates to the 
impact parameter. The results are very similar to theory predictions in case of 
crossover from hadron gas to quark gluon plasma. According to our analysis, 
the possible crossover should start at ~ 6 and end at ~ 24. 
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i 1. Introduction 



2 Some of the most important questions about the transition to the quark 

3 gluon plasma (QGP), a new state of matter with partonic degrees freedom, 

4 are not yet fully answered. Among them the location of phase boundaries 

5 between hadronic gas and the QGP. The results of lattice QCD simulations 
e concerning the order of phase transition depend strongly on the number 

7 of quark flavors and on the quark masses. For vanishing baryon chemical 

8 potential //& = 0, the nature of transition can be a genuine phase transition 

9 (first order or continuous), or just a rapid change (crossover) over a small 

10 temperature range jl|. Estimates of energy densities which can be achieved 

11 in ultra-relativistic pp or pp collisions with high multiplicities suggest values 

12 sufficiently high for experimental formation of the QGP {2]. 

13 However it may be that, unlike what happens in heavy ion interactions, in 

14 pp and pp the central blob of created matter never thermalizes yl, although 

nnn 

15 there are different opinions |2J, |4H7j which predict that thermodynamics con- 

16 cepts may be applied in pp or pp high multiplicity events. 

17 Probes of equation of state are among possible signatures of phase tran- 
ib sition or crossover. The basic idea behind this class of signatures is the 

19 identification of modifications in the dependence of energy density e, pres- 

20 sure P and entropy density a of hadronic matter on temperature T. One 

21 wants to search for a rapid rise in the effective number of degrees of freedom, 

22 as expressed by the ratio e/T 4 or a/T 3 , over a small temperature range. One 

23 can expect a step-like rise as predicted by lattice simulations (Fig. [I]), more or 

24 less steep depending from the presence of transition or crossover, and from 

25 the order of the transition in the former case. Finite volume effects may 
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26 cause important consequences for e/T 4 and cr/T 3 : the latent heat and the 



27 jump in the entropy density are considerably reduced for small systems 



8]. 



28 Besides that, the critical temperature may shift to higher temperatures and 

29 the width of the transition may broaden for smaller volumes and there may 

30 be a smoothening of singularities due to the finite size of the system js-lO]. 



31 In 1982 it has been suggested by Van Hove [12| that an anomalous be- 

32 havior of average transverse momentum (pt) as function of the multiplicity 

33 could be a signal for the occurrence of a phase transition in hadronic matter. 

34 His conjecture is based on the idea that the (pt) distribution of secondaries 

35 reflects the temperature of the system and its evolution in the transverse 

36 direction, while the multiplicity per unit rapidity provides a measure of en- 



tropy 



14|. In a recent paper [ij]] one of us showed that from 22 to 



38 7000 GeV in 21 (p?) vs pseudorapidity density curves there is a slope 

39 change at = 5.5 ± 1.2. Signals related to these slope changes may indi- 

40 cate transition to a new mechanism of particle production. Many years ago, 



4i m 



161 ]. we pointed out that pp at ISR and pp data at CERN collider showed 



a kind of jump at = 6 and that it had to be investigated as a possible 



dij 

phase transition signal [17|. In 2002, Alexopoulos et al. 



1* 



assumed that 



44 the system produced in pp at v/i = 1800 GeV for ^ > 6.75 was above the 

45 deconfinement transition to explain their experimental results. 

46 In present article, taking into account experimental results in pp and pp 

47 at high energies 19l432l|. we show how measured physical quantities satisfy 

48 relations which, given proper approximations and correspondences, can give 

49 a representation of the equations of state (EOS) that describe the created 

50 system in the central region in pseudorapidity in high energy pp and pp 
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Figure 1: Results of the lattice simulations of QCD for T > T c (critical tem- 
perature) and from ideal hadron-gas model for T <T C . 
[Taj entropy density o scaled by T 3 calculated in the hadron-gas model and 
by lattice simulations of QCD shown as function of temperature. The verti- 
cal line indicates the critical temperature. 

[Tbl sound velocity c^ shown as function of the energy density e. 

ITcl temperature dependence of the square of the sound velocity at zero 

baryon density as function of T. In this case the critical temperature T 

is equal to 170 MeV. 

From 11 1 . ^ 



51 collisions. Starting from (p T ) vs experimental results together with the 

52 estimation of the size S of the interaction area, which is obtained from the 

53 measurements of the radii of emission in function of multiplicities 23l-|26| . 

54 or from a model which relates multiplicity to impact parameter we 

55 obtain relations among (pt) and particle density as, which seem to resemble 

56 EOS curves predicted for hadronic matter with crossover to QGP. The (pt) 

57 and S vs relations contain the relevant information, which translates in 
5B (p T ) vs as correlations. 

59 According to our knowledge, this is the first attempt to obtain an esti- 

eo mation of the complete EOS for hadronic matter using experimental data 

ei only. 

62 2. Methods 

63 The experimental results and the approximations made in this work are 

64 the following. 

65 2.1. (p T ) VS^ 

66 As we mentioned before, (pr) vs correlation at about = 6 shows 

67 a slope change in all the experiments. In Van Hove scheme (pt) reflects fem- 
es perature and the system evolution. On the other hand, the biggest part of 

69 emitted particles is constituted by pions and the pion (p?) is rather insensitive 

70 to flow [34]. Thus, not identified charged particles (px) may be considered 

71 as an estimation of the system temperature because it's not influenced very 

72 much by a possible transverse expansion. Furthermore, transverse radius 

73 Rside 

vs pair transverse momentum in pp Bose Einstein correlation mea- 



74 sures 



24J show that, at least until ~ 3.4, results are consistent with the 
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absence of transverse expansion, which further supports the adoption of (p?) 
as an identifier of temperature, because it's little affected by the expansion. 

Since a substantial number of pions is the product of resonance decay 
and the particles originating from the resonance decays populate the low px 



79 region 



11], in this work we consider mainly (px) vs correlations with 



so a pt min cut (> 400 MeV/c in CDF experiments Run I and Run II , > 500 

si MeV/c in ALICE, and two different cuts, > 500 MeV/c and > 2500 MeV/c, 

82 in ATLAS experiment) in order to work with (p?) values less influenced 

83 by this effect. Furthermore, diffractive events are substantially reduced for 



84 ^ > 2 in (p T ) vs ^ plots with p T min > 400 MeV/c [2l|, |22|. 

85 We will show anyway also some results for (px) computed with pt mm cut 
se and > 100 MeV/c. The structure of the relations we are going to show is 
87 still present in these measures. 

sa In this work, (pj) computed for different pt cuts will be plotted without 

89 the application of corrections due to the cut in the used pt range, apart 

90 from the case of events energy density estimation, in which we will use a 

91 corrected (pt)- Regarding —^ L , it is computed from the number of particle 

92 in a given region of pseudorapidity r\ and px, dividing by the amplitude of 

93 the 7] range and properly correcting for p? cuts. In order to perform this last 

94 correction, we considered curves for the different experiments, measured 

95 with and without px cuts, and multiplied by the ratio between correspondent 

96 values of in the central region. All data are obtained from minimum 

97 bias experiments. For CDF run II 1960 GeV, high multiplicity trigger data 

98 are added to minimum bias data, for charged particle multiplicity N r h > 22 



( \t]\ < 1, Pt > 400 MeV/c corresponding to corrected > 22 
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100 Where available, we considered also raw data results (i.e. computed with- 

101 out experimental inefficiencies corrections) because the (pr) vs plot and 

102 its derived plots are much sensible to experimental losses and, on the other 

103 hand, the application of corrections may involve some "smearing" of data 

104 which could highly modify the analyzed effects 19|, [35|. For reasons of space 



105 we don't show behaviors for raw data in this paper, because results are very 
loe similar to those for corrected data. 

107 2.2. Entropy Density Estimation 

The initial energy density in the rest system of a head-on collision has 
been argued to be [6(: 

Hit) 



dr] 



V 

V denotes the volume into which the energy is deposited. Similarly the initial 
is Q: 



entropy density 

dN ch _ 3 
dn 2 

a ~ 



V 

los As a result, e is equal to a ■ (pr)- The volume V may be estimated as 
ioo V = S ■ ct, where S is the interaction area and ct is a longitudinal dimension 
no we can traditionally consider to be about 1 fm long. 

In order to study our system, we will use the quantity 

dN ch _ 3 
dn 2 

in as an estimation of entropy density. In models like color glass condensate and 
ii2 percolation, the system physics depends on cr s |36l-K 



For the estimation of 



n3 the area of interaction S, we proceed in different ways, our target being the 
H4 obtainment of results which are robust respect to the definition of the area. 



115 On the other hand, we are more interested in relations between variables 
lie than in their absolute values. 



2.3. Bose Einstein correlation for emission region size estimation 

Using Bose Einstein correlation among emitted particles, measurements 



ii9 of p article emission regions in many pp and pp experiments have been done 



26 



31 



32 



39 



23 



4l| In [23|, |24|, as already mentioned, the measurement of R S id e 



121 in function of in pp, shows that the transverse radius doesn't depend on 

122 krp for low ^ values ( < 3.4). This can be explained by the absence of ex- 

123 pansion of the particle emission source, at least at these values. For ^ sk - 

124 values greater than 7, there is a dependence of R s ;d e on fc^, so that probably 

125 a source expansion is possible at least from this value. It is thus possi- 

126 ble that a new phenomenon is started in events with between 3.4 and 

127 7. The hypothesis of no expansion for low values lets us approximate 

128 the initial interaction section radius to be coincident to the final emission 

129 radius. We take into account that for values greater than about 7.5, 

130 this approximation is more uncertain. Furthermore, resonance effects are 

131 present, but, once more, we are not interested in the absolute values of the 

132 interaction section, but in its behavior in function of Not taking into 

133 account these effects yields a systematic error on the value of the radius, 

134 which we consider invariant for different values of ^ iL - Given the similarity 

135 in both the behavior and the absolute values of R S ide and invariant radius 

136 R inv versus multiplicity, we use R inv as an estimation of the interaction region 

137 radius, mainly because R nv data where measured for a larger range than 



138 R S ide ones 
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-26 



31 



32 



39 



40j. In Fig. R inv is shown CIS db function of 



pseudorapidity density. In the left we only show data for CMS (preliminary) 



140 and ALICE (preliminary) at (k T ) ~ 0.35, while in the right side we show 
mi the same results along data from other experiments (UAl, ABCDHW ISR, 



142 STAR). We fitted the data of Fig. [2^ with two functional re 



ations between 



143 R inv and the first is linear in the cube root of [23|, Ul, |42[ and the 

144 second is linear in cube root of for > 7.5, matched with a 5 th de- 

ar; dr] 1 

145 gree polynomial fit for smaller values. The first fit gives a N p oF = 0.84 

2 

we with p-value: 0.67, while the second gives a N p F = 0.45 with p-value: 0.95. 

147 Considering these results, we opted to use the second fit for the following 

us analysis. 

149 It has been stated that the behavior of radii in function of doesn't 

i — ~\ c "? 

150 depend on the experiment energy [41]. Data in Fig. [2] seem to confirm this 

151 statement, and justify our choice of a single relation for R inv vs for all 

152 energies. 

153 In order to estimate the interaction region, we used the following alter- 

154 natives: 

155 1. An area obtained using R; nv from the combined (polynomial + linear 

156 in cube root of —^) fit from Fig|2al This choice may overestimate the 

157 interaction region in case of system expansion, being R inv a measure of 

158 the emission region; 

wo 2. following ALICE results in R S id e vs hp, we make the hypothesis that 
wo no expansion is present in events with sufficiently low -^ sk - So we 

lei use Ri nv from the left (polynomial) part of the combined fit in Fig [2^, 

162 then we use a constant radius for > 7.5 as an estimation of the 

dr] 

163 dimensions of the initial region before the possible expansion, making 

164 the assumption that at ~ 7.5 the interaction region reaches its 
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23 



26|, 



31 



32 



40|. Left 



plot shows data from CMS (preliminary) and ALICE (preliminary) both at 



(kx) — 0.35, along with linear fit in 



dN ch 
drj 



and a combined fit: 5 degree 



polynomial for < 7.5 matched to linear fit in d x 



h 3 



for 



dN r , 



> 7.5. Right 



drj ' drj drj 

plot shows data from the first plot along with data from UAl, ABCDHW 
ISR, and STAR experiments. 
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165 maximum; at = 7.5 the Ri nv value is 1.08 fm; 

lee 3. an area obtained from a model which relates the impact parameter to 

167 the multiplicity of events [331 ] . 

168 From values and from interaction areas, estimated as described 
leg above, we obtained the values of density of particles for transverse area. 

170 We considered the simplified case where the central blob volume V is the 

ni same in all collisions for a given [12]. We estimate an average a from 

172 the ratio between and the estimated average V. 

173 2.4. a s / (pt) 3 vs (p T ) 

174 Using the estimated 05, the relation (p T ) vs as can be studied. A slope 

175 change in (p?) vs as plots is found at as between 2.5 and 3 fm~ 2 , depending 

176 on the method used for the estimation of area S and corresponds directly to 

177 the slope change seen in (p T ) vs at ~ 6. 

178 Starting from as and (pr), we plotted as/(pr) 3 vs (pt) curves, as an 

179 experimental approximation of a/T 3 vs T curves. See Figs. [3] and HI We 

180 obtained v ery similar as/(pr) 3 vs (pt) curves from other pp and pp experi- 



i8i ments 



27 



29|-|3jJ (not shown). 



182 In figures, we put labels with corresponding values for interesting 

183 points, in order to relate these points to the characteristic values in —f^- 

184 In the different plots, different regions are recognizable. In particular, in all 

185 plots we see that from the as/{pr) 3 value corresponding to ~ 2, up to 



chj 

186 a value correspondent to ~ 6, the curve is almost flat, then rises very 



187 quickly. This behavior is similar to the one in a/T 3 curve, in presence of 

188 crossover, starting from a state of matter, identified by a/T 3 nearly constant 
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Figure 3: cts/{Pt) vs (p T ) plots. Area S from "5 degree polynomial + 
constant after ^ > 7.5" fit. 

ALICE at yfl = 0.9 TeV, 0.5 < p T < 4 GeV/c, \r]\ < 0.8, Minimum Bias. 
ATLAS at yfe = 0.9 TeV, p T > 0.5 GeV/c, \r}\ < 2.5, Minimum Bias. 
3 ATLAS at = 7 TeV, p T > 0.5 GeV/c, M < 2.5, Minimum Bias. 
ATLAS at y/Ts = 7 TeV, p T > 2.5 GeV/c, M < 2.5, Minimum Bias. 
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Figure 4: a,s/ (pt) 3 vs (pt) plots. Area S from "5 th degree polynomial + 



constant after ^ > 7.5" fit. 

dr] 



1 CDF Run I at a/s = 1.8 TeV, p T > 0.4 GeV/c, < 1.0, Minimum Bias. 

CDF Run II at v/i = 1.96 TeV, p T > 0.4 GeV/c, \rj\ < 1.0, Minimum 
Bias + High multiplicity trigger. 

3 ATLAS at yfe = 7 TeV, p T > 0.1 GeV/c, < 2.5, Minimum Bias. 
CMS at v/i = 7 TeV, p T > GeV/c, \r}\ < 2.4. Minimum Bias. 
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189 (region 2 < < 6), and a crossover starting at ~ 6 (Fig. CP . Besides, 

wo in plots with many points at high ^ values (ATLAS with p T > 2500 

191 MeV/c, and CDF Run II 1960 GeV with p T > 400 MeV), we observe a 

192 strong slope change around corresponding values of about 24 or higher. 

193 It's worth noting that what seems to be a different behavior in the left side 

194 for ATLAS with p T > 2500 MeV/c (Fig. [33]), is only due to the fact that all 

195 points correspond to > 7, apart from the first point, which correspond 
we to ~ 3.4. The ratio between a si (pt) 3 values corresponding to > 24 

197 and those corresponding to < 6 varies from 2 to 3, depending on the 

198 area calculation method used for the estimation of as- This ratio in the case 

199 of EOS would correspond to the ratio between the number of the degrees of 

200 freedom of the state before and after the transition or the crossover. We note 



that for small size systems as it wou 



d be in the pp case, the jump in entropy 



density is considerably reduced 8Nl0j in comparison to the theoretical infinite 



volume case. In plots with p T > 100 MeV/c (ATLAS 7 TeV) or p T > (CMS 
7 TeV), the first points have constant (px) with varying which leads 

to an initial steep rise. After that, the curves assume the same behavior of 
previously seen plots. 



2.5. Sound velocity c 



209 



One of the physical quantities used to characterize the state of a system is 

(j dT 

its squared sound velocity, defined as c^ = — • — — , for constant V [11]. In our 

T da 

210 study, we approximate it with c^ = S ■ ^ f T . It is really interesting that 

(Pv) da s 

2n if (pt) is proportional to T and if as is proportional to the entropy density, 
212 then the Cg value obtained in this approximation is equal to the right value 
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213 of Cg = ^- ■ — — , because proportionality constants cancel out. In order to 

T da 

214 obtain our c 2 estimation, from (p?) vs curves and from as values, we 

215 compute the curve (pt) vs as, to which we apply numerical derivation. We 

216 cope with the statistical fluctuation in data points using a combination of 

217 Gaussian and Savitzky-Golay filters 43|. Examples of cf vs (pr) curves are 

218 shown in Figs. EH and EcJ 

219 The so obtained c 2 estimation resembles the typical shape of a phase 

220 transition or a crossover: a descent, a minimum region and a following rise, 

221 as it's also obtained analytically from EOSs which present a phase transition 

222 or a crossover. The minimum value reached by the estimation of c;? in the 

223 different experimental curves varies from 0J)8 to 0.18 and could correspond 

224 to what it's called the 

225 Recently Refs. 



44 



)OS softest point [17] . 



46] estimate c 2 minimum value for realistic EOS to 



226 be around 0.14. From eg ~ (pt) ■ as we compute c 2 vs es curves, that are 

227 approximations of c 2 vs energy density. We report these curves in Figs. [5b] 

228 and l5dl 

229 In this case, values are calculated using (p T ) values from (p T ) vs 

230 curves with no px m in cut at corresponding energies, estimated in correspon- 

231 dence with the different values. As Figs. [5] show, the numerical estima- 

232 tion of Cg vs es, is characterized by a maximum at low energy density fol- 

233 lowed by a minimum region, which is obtained for values in range 1.5 — 2.0 

234 GeV/fm 2 , and a subsequent rise. 

235 We note that as computed here is an estimation of the energy density 

236 for pseudorapidity unit and unit of transverse area. In order to estimate the 

237 volume energy density this should be divided by ct. 
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Figure 5: c s 2 = — — ■ — 

\Pt) da s 



(d) 



vs (pt) or es, using two different fits for area S 



estimation. 

ALICE at v/i = 0.9 TeV, 0.5 < p T < 4 GeV/c, \rj\ < 0.8, Minimum Bias; 
area S from "5 th degree polynomial + linear in after > 7.5" fit. 



a) Cg vs (j>r),[5bl vs e^. 
CDF Run II at y/s = 1.96 TeV, p T > 0.4 GeV/c, |r?| < 1.0, Minimum Bias 
+ High multiplicity trigger; area S from "5 th degree polynomial + constant 



after ^ > 7.5" fit. 

drj 



3 c 2 s vs (p T ), 



c 2 s vs eg. 
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238 3. Discussion 

239 The shape of the <Js/{pt) 3 approximation to the EOS is very similar, 

240 using both R inv from the fit on all space and R inv fitted up to = 

241 7.5 and then maintained constant. It slightly varies when using the area 

242 from the impact parameter model, but the slope change at <Js/(pt) 3 values 

243 corresponding to around 6 is still present, as well as the change at 

244 cts/(pt) 3 values corresponding to about 24. 

245 In order to avoid possible systematics due to calculation involved in the 

246 area definition, we plotted directly ^V(pr) 3 vs (pt)'- this is equivalent to 

247 obtain erg/ (pt) 3 curves considering a transverse section which is constant for 

248 all values. For space reason we don't show these plots in this paper. 

249 In this case the shape doesn't resemble an EOS shape anymore, but the 

250 slope changes at ~ 6 and ~ 24 are still present, because they are 

251 contained in the (p?) vs correlation. 

252 The shape of the curves obtained from experimental data (as / (pt) 3 vs (pt) , 

253 Cg vs (pt) and Cg vs energy) depends on experimental (pr) vs curves and 

254 from the value of the area used to obtain density sigmas. Systematic errors 

255 in (pr), an d Rinv measurements don't lead to appreciable variations in 

256 (pt) vs as behavior, which is what we are interested on. 

257 It seems to us that the main result of this work is that putting together 

258 experimental data of (pr) vs and Ri nv vs curves are obtained which 

259 are the reproduction of theoretical EOS curves. 

260 Regarding model comparison, we obtained as/ (pr) 3 vs (pr) plots starting 

261 from Montecarlo curves (Pythia ATLAS AMBT1 and Pythia8 for ATLAS 

262 and CMS experiment respectively), which are shown in Fig. [6] 
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Figure 6: cts/(pt) 3 vs (pt)- Comparison with models. 

6aJ ATLAS at ^ = 0.9 TeV, p T > 0.5 GeV/c, \rj\ < 2.5, Minimum Bias and 
Pythia ATLAS AMBT1. 

6b] ATLAS at v/i = 7 TeV, p T > 0.5 GeV/c, H < 2.5, Minimum Bias and 
Pythia ATLAS AMBT1. 

M ATLAS at y/s = 7 TeV, p T > 0.1 GeV/c, |^| < 2.5, Minimum Bias and 
Pythia ATLAS AMBT1. 

6dJ CMS at v/i = 7 TeV, p T > GeV/c, < 2.4, Minimum Bias and 
Pythia8. 

Bel ATLAS at yfi = 7 TeV, p T > 2.5 GeV/c, \r]\ < 2.5, Minimum Bias and 
Pythia ATLAS AMBT1. 
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263 Some models on which tuning has been done, for example with CDF Run 

264 II data at 1960 GeV for p T > 400 MeV/c, well reproduce the (p T ) vs 

265 curve at higher (7 TeV) or lower (0.9 TeV) energies with pt > 500MeV/c. It 

266 is clear that in these cases, starting from the (p?) vs curves of models and 

267 using Bose Einstein correlation or the impact parameter-multiplicity relation 

268 for as estimation, curves similar to the experimental ones are obtained. On 

269 the other hand, models don't predict well (px) vs curves with low pr mm, 

270 and consequently <Js/(Pt) 3 vs (pt) curves as shown for the comparison of 

271 models at 7 TeV for CMS and ATLAS data, respectively with pt > MeV/c 

272 and p T > 100 MeV/c 20 J3, 3- 

273 The interpretation of curve shapes as experimental "estimation" of EOS 

274 depends on how much likely are the correspondences between (pt) and T, 



275 and between measured as and entropy. 



276 4. Conclusion 



The result we consider to be the most important is the following: in 



278 many experiments 
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27 



3 if ] from 31 GeV to 7000 GeV, starting from 



279 (pt) vs and using results from measures of radii with Bose Einstein 

280 correlation or from a model that relates impact parameter and multiplicity, 

281 we obtained that a s / (pr) 3 vs (pt) and Cg = 5 ■ reproduce the shape 

(p T ) da s 

282 of hadronic matter EOSs and squared sound velocity respectively, in presence 

283 of crossover or phase transition. From the plots, a neat change around dA "' 



dr] 

around 6, where the crossover or the phase transition seems to start, and 
another possible change at around 24 are observed. The curve c^ vs es 
has a minimum around a "transverse" energy density of about 1.5 GeV/fm 2 . 
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28? In order to understand if these behaviors have a real physical meaning 

288 or are just casual, results of measures in the following regions should be 

289 compared: 2 < ^ < 6, ^ > 6, 6 < ^ < 24 and ^ > 24. 
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